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ABSTRACT 

A semi-empirical control-volume approach to ejector modeling for transient performance 
prediction is presented. This new approach is motivated by the need for a predictive real-time 
ejector sub-system simulation for STOVL integrated flight and propulsion controls design 
applications. Emphasis is placed on discussion of the approximate characterization of the 
mixing process central to thrust augmenting ejector operation. The proposed ejector model 
suggests transient flow predictions are possible with a model based on steady-flow data. A 
practical test case is presented to illustrate model calibration. 

INTRODUCTION 


For the short take-off vertical landing (STOVL) aircraft configuration shown in Figure 1, the 
thrust augmenting ejectors placed at each wing root must be included in the propulsion system 
simulation. Integrated flight and propulsion control (IFPC) studies require ejector component 
simulations that must run accurately in "real-time". Since accuracy is typically synonymous 
with a high level of model detail, constructing a transient ejector model to fit the simulation 
requirement is not an easy task. 

A useful ejector characteristic to know the frequency response of ejector thrust. If the ejector 
response is outside the bandwidth of the flight control system then a simple quasi-steady 
ejector model or look-up table is all that is needed. On the other hand, if the ejector response 
falls within the range of the flight control (or if a new ejector design has an unknown response) 
a high-fidelity ejector simulation is required; this paper explores an ejector modeling approach 
to fill the latter requirement. 


There are two basic perspectives generally employed to set the mathematical framework for 
the ejector mass, momentum, and energy balances; the choice is easy when the intended 
simulation application is considered. One approach is based on first principles and typically 
blends the full unsteady Navier-Stokes equation with a tubulence model. The result is cast 
in, for example, finite-difference form. Although this approach provides a desireable ejector 
design flavor to the analysis, the solution to the final set or equations is not amicable to real-time 
execution in the foreseeable future; another problem concerns identification of the appro- 
priate turbulence model to use for description of the mixing region. 
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A second, more analysis oriented approach, involves a control-volume formulation. Highly 
detailed descriptions of the mixing region shear layer can be exchanged for a semi-empirical 
model and the potential to run in real-time. In this paper, the intended ejector application as 
a force effector leads to tailoring of the semi-empirical model to match thrust predictions 
(matching is not primarily for the field variables). This compromise is acceptable in light of 
the simplicity the model offers. 

An attractive feature of a steady-flow control-volume ejector analysis is the ability to collapse 
the complex primary and secondary flow mixing process for steady flow into one unknown 
quantity - the mixing effectiveness parameter (3 . Although ideal mixing yields a mixing 
parameter of unity, the more typical non-uniform flow profile (non-ideal mixing combined 
with boundary-layer effects) in real ejector configurations generally produces mixing 
parameter values between 1.05 and 1.2. The key in a useful control-volume analysis is knowing 
what the exact mixing effectiveness value is for a given ejector configuration. In the analysis 
to follow it is assumed steady-state ejector data is available so that [3 can be backed out. 

Under the premise the steady-state ejector model has been "calibrated" (that is, [3 and all 
surface integrals are known), the time-dependent volume-integrals are the only remaining 
unknowns to be treated in the control volume analysis of an unsteady flow. In the present 
work, sub-dividing the mixing region (streamwise) into a finite number of subvolumes allows 
Abramovich-type self-similar turbulent flow profiles to conveniently link the local primary 
and secondary field variables. The result is a set of equations for each mixing region sub- 
element containing the desired field variable time derivatives. From this analytic springboard 
a kinetic energy exchange function is introduced to complete the interaction between the 
streams and provide closure for the mixing region analysis. 

MIXING REGION REPRESENTATION 

Numerous finite-volume descriptions for steady-state ejector flows are available in the lit- 
erature (Drummond[4], Porter and Squires[6], Addy and Dutton[2]). For simplicity in analysis, 
ejectors are typically divided into inlet, diffuser and mixing regions. Since inlet and diffuser 
dynamics can, in the first approximation, be assumed away with quasi-steady approximations, 
only the ejector mixing region shown in Fig. 2 is of interest in the present work. A control- 
volume mixing region analysis generally involves the one-dimensional mass, momentum and 
energy equations below: 
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These equations reflect the use of the usual assumptions of adiabatic frictionless flow (at the 
ejector wall) [4]. A steady-flow analysis involves only the surface integrals on the right- 
hand-side of Eqns.( l)-(4). It is the time derivatives of the volume integrals on the left-hand-side 
of each equation that characterize the unknown dynamics of the ejector mixing region. 

Remark on the Energy Equation 

Equation (3) represents the energy equation derived from the first law of thermodynamics -- 
this form of the energy balance is known as the general energy equation. The general energy 
equation is a balance of energy due to heat and work. The product of the local fluid velocity 
and the momentum equation (Eqn.2) yields the mechanical energy equation. Subtracting the 
mechanical energy equation from the general energy equation produces the heat equation , 
the latter of which is a simple statement of conservation of heat content in the ejector mixing 
region, 


JtiA h 'f) dV ' -£ phv ' ncM 


It is important to note in the present work that the general energy equation is split into its 
mechanical and thermal energy components. Although Eqn.(4) is retained to represent the 
heat balance, a semi-empirical kinetic energy exchange model is substituted for the mechanical 
energy equation (the details of the function are discussed later). The idea to emphasize here 
is that an emprically-based transfer function is proposed to be a satisfactory replacement for 
the detailed mechanical energy modeling that would otherwise be required for the turbulent 
flow characterization. 


Field Variable Approximation 

Extensive experimental data supports the basic Abramovich[l] non-dimensional field variable 
representations of the co-flowing primary and secondary streams. General features of the 
ejector jet geometric approximation are shown in Fig. 3. Extending from the mixing region 
inlet plane there exists a potential-core region characterized by a fairly uniform centerline 
velocity, with no transverse component. This is distinguished from the mixed- flow region where 
the centerline velocity decay arises from momentum transport to the entrained fluid. In the 
characterization of an element of the ejector mixing region, it is assumed the velocity can be 
approximated by the following 2-D planar turbulent jet self-similar profiles for co-flowing jets: 


where 


v = 




i/.a-*)-!/.,* = /(<*>) 

( 5 ) 

_ /(l-^- 5 ) 2 ,0<^<l\ 

\ o, i<^<i / 

(6) 

x/b, | = B/b 

(7) 


3 


Inclusion of the core region requires a slight modification of these profiles; details are given 
in Abramovichfl]. Since there is a static pressure matching condition at the jet boundary in 
the ejector mixing region, a uniform transverse pressure distibution must be assumed; in the 
longitudinal (axial) direction, however, a finite pressure gradient exists - this feature partially 
distinguishes free-jet and confined-jet analyses. Although the general ejector model admits a 
function A to characterize the normalized pressure distribution, the uniform transverse 
pressure suggests 


-Z—Z — = A($) = [1.0<S<|] (8) 

Pm P e 

To reduce confusion in the use of the ideal gas law, pressure and density should have similar 
profiles, 


P = P.(l-A) + Pm A 


(9) 


Since these functions imply a uniform temperature profile, a compromise has been introduced 
in the mixing region representation for temperature; Abramovitch suggests the non- 
dimensional temperature profile be approximated by the square root of the velocity profile. 

Finite Volume Descretization 

Figure 4 illustrates the finite volume descretization employed in the present work. Since the 
implicit assumption has been the self-similar profiles link the primary streams in the transverse 
direction, descretization occurs only in the streamwise direction. Because the characteristic 
streamwise velocity is generally "high", very few volumes are needed to capture basic flow 
effects. Computations in the present work were performed with the mixing region divided 
into 5 elements. 

The combined effect of the analytic interstream relationships and the streamwise descreti- 
zation is to provide a virtual computational grid, also shown in Fig.4. 

Application to the Conservation Equations 

Only a summary of the application of the non-dimensional field variable profiles to the control 
volume conservation equations is presented here since details can be found in reference [4]. 


Conservation of Mass 

For the generic sub-region k, bounded by surfaces at * and j, mass conservation yields 
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First, the mass flux is written 
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From which integration over the self-similar profiles yields 


m, = 2W'b lP< (0.45u in + 0.55t/, + (i- l)u e ), = 2Wb i Z l (12) 


The time-derivative term in the mass conservation equation is given by 
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where the characteristic density for the finite volume is now approximated by the value of the 
density at station j (this is a reasonable assumption as long as field variable gradients are 
"modest" in size). If the characteristic jet expansion width, b, also assumes its value at j, then 
substitution and re-arrangement of the continuity equation yields 
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Computation of the jet half- width, b, derives from an incompressible flow assumption applied 
to the first finite volume; a rectilinear jet expansion is assumed therefrom and compressible 
flow restored for subsequent calculations. 

Momentum Equation 

The procedure described for the continuity equation can be applied to the momentum 
equation. Omitting the details for brevity, the result for the centerline velocity derivitive is 
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Energy Equation 

Introduction of an ideal gas assumption simplifies the heat equation to the form 


d 

dt . 
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Non-dimensional velocity and pressure profiles provide the pressure derivative result 
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where Zi has been derived previously and Zq is given at each station by 
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Equations (14), (15), and (18) are the three differential equations required for computing the 
time derivatives of gas density, pressure, and jet centerline velocity in the ejector mixing region. 
Integration of these equations can be combined with a quasi-steady diffuser analysis to predict 
the ejector discharge conditions as a function of time. Ejector thrust is approximated by 


T = 
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where p is the mixing effectiveness parameter discussed in the introduction. 

Closure of the initial value problem for the mixing region requires an analytic approximation 
for the mechanical energy removed in the formulation of Eqn.(4). The purpose of the next 
section is to propose an appropriate function. 

Entrained Kinetic Energy Approximation 

This section provides an approximation for the turbulent flow kinetic energy exchange 
mechanism to characterize the influence of primary flow changes on the secondary flow. A 
discussion of the proposed kinetic energy balance is followed by clarification of the primary 
and secondary kinetic energy representations. 

Kinetic Energy Balance 

Computations for a specified steady-state condition show that the change in kinetic energy 
due to mixing is not the same for the secondary flow as it is for the primary. In fact, the gain 
in kinetic energy of the secondary flow is entirely due to the mixing process, while the mixing 
loss of the primary flow is only a fraction of its total loss. In balance, however, the total change 
of kinetic energy of the primary flow is greater than that of the secondary flow. 

In the work of Korst and Chow[5] the relationship between the change in entrained flow 
kinetic energy and the total primary flow kinetic energy for a shear layer is given in functional 
form by 


A KE XS = f[kE xp ,^ (21) 

where a reasonable match between theory and experiment is given by 

o = 12(1+0.23 M lP ) (22) 
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The difficulty with the kinetic energy function as given above is that it represents a quasi-steady 
constant-pressure flow approximation and therefore cannot be used in its present form for 
the transient flow analysis. To entertain local transport of energy between the primary and 
secondary flows, consider the change in secondary flow to be a combination of changes in 
primary and secondary flow kinetic energies due to mixing alone, 

AKE ls = G^AKE lPtm ,AKE XStm ,^j (23) 

where the subscript m denotes the change in kinetic energy due exclusively to mixing. In 
conjunction with numerical experiments, a specific form for the function has been determined 
to be 


A KE 


t + At 
IS 


AKE 


1 S,m 


+ A K E [ P m C 



(24) 


Note the introduction of this engineering approximation also results in the introduction of an 
undetermined constant, C 1# This value is established in the present work by matching transient 
solution asymptotes with steady-state performance data; an example is presented later a the 
case study. The alternative is to establish N computations of the kinetic energy exchange to 
coincide with the N control volumes of the mixing region; the present method is more rapid 
since it permits post-processing of information at the completion of mixing region calculations. 

Knowledge of the total change in secondary flow kinetic energy permits updates to the sec- 
ondary flow state during integration of the field variable time derivatives. After the velocity 
of the secondary flow has been extracted from the new secondary flow kinetic energy, the new 
inlet density and pressure are computed from an ideal flow inlet analysis. Brief descriptions 
of the mixing kinetic energies are given below. 


Secondary Stream Energy Change 

Computation of the gain in secondary flow kinetic energy due to mixing is given by [5] 


A KE = 



(25) 


where <; defines the jet boundary streamline illustrated in Fig.5 (the streamline is positioned 
so that the secondary mass flow through station i is equal to secondary mass flow through j). 
For the present discussion the dividing streamline position, b , is assumed known; Korst and 
Chow[5] discuss the typical approach of analysis. Expanding the equation for the change in 
kinetic energy yields 


A KE = 1/6 p 


(v 3 -vv* 



(t > 3 -vv z e )d% 


(26) 


Substitution of the self-similar profiles into this expression and integrating the result provides 
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where the detailed expressions for the integrals Hj and F[ are given by Drummond[4]. What 
is important about the integrals is that they are independent of time and have closed form 
algebraic solutions; they also need only be computed once for a given shear boundary layer 
value. 


Primary Flow Energy Change 

Similar to the way in which the change in secondary flow kinetic energy was computed, the 
energy loss of the primary flow is given by 


A KE 



(28) 


where the limits of integration reflect interest in the domain of the primary jet cross-section. 
Completing the integral for station i yields 
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Again, reference [4] provides the expressions for the integrals implied in F( and /// 

DISCUSSION OF CASE STUDY 

For complete verification of the proposed ejector modeling technique it is necessary to have 
available transient ejector performance data to compare with calibrated simulation output. 
At the present time, however, only steady-state tests have been conducted at NASA Lewis 
for the STOVL ejector shown in Fig.l (see summary of Corsiglia et. al.[3]), but this is proposed 
to be adequate for ejector model calibration; by "calibration" it is meant that the unknown 
mixing parameter f3 and the kinetic energy constant C ! have been computed for a given 
ejector configuration. In this regard a practical STOVL ejector test case to examine involves 
a step-function in the primary nozzle flowrate. 

The first constant to compute is the mixing effectiveness parameter, (3 . Here, |3 is selected 
as the value that appears to provide reasonably accurate thrust predictions for all steady flow 
data points within a given performance window. Tests of the E7-D ejector (Fig.l) covering 
subsonic and choked primary flows suggested (3 = 1.055 provides the best general match 
between theory and experiment. 

For a step-function change in the primary flow, the ejector, initially operating at a steady-state 
condition, asymptotically approaches a second steady-state condition. The desired Cj value 
is the one which provides a match between the predicted asymptote and the known thrust at 
the second state. Figure 6 illustrates the situation for subsonic primary nozzle flow (within 
the performance window for (3 above) experiencing a change to 48.17 kg/s from 41.23 kg/s; 
the corresponding thrust levels are 212 N and 270 N, respectively. Thrust predictions for 
various Cj values are shown in Fig. 6; the "best" asymptotic match appears to be for Cj between 
0.3 and 0.35. 
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The thrust predictions shown in Fig.6 appear to have a second-order response to the primary 
flow step-function. It seems the level of C j primarily influences the level of the thrust rather 
than shift the peak thrust overshoot. Under the second-order assumption the damping ratio 
is approximately 0.75 and the natural frequency of the ejector is on the order of 300 Hz. 

CONCLUSION 

A mathematical ejector model that includes system transients and has the potential to run in 
real time has been presented. The finite volume method permits time derivative of the mixing 
region field variables to be rapidly evaluated. Although final code verification awaits the 
availability of transient ejector data, ejector thrust predictions are intuitively reasonable. 
Future work must explore the proposed kinetic energy transfer function within a larger ejector 
performance window. 
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NOMENCLATURE 


A cross-sectional area 
b jet half-width 
B channel width 

C j coefficient defined by equation(24) 

KE kinetic energy 

rii mass flowrate 

m cv control volume mass 

M momentum flux 

M cv control volume momentum 

M i fMach number of primary nozzle flow 

P pressure 

t time 

v velocity 

x transverse coordinate 
streamwise coordinate 


Z i integral of defined function 
(3 mixing effectiveness parameter 
Y specific heat ratio 
A normalized pressure profile 
p density 

4> function defined by equation(4) 
a turbulent flow constant, equation(20) 

^ dimensionless transverse coordinate 

t; dimensionless jet streamline boundary 
subscripts 

cv control volume 
g entrained stream 
m primary stream centerline 
1 s secondary steam at station 1 
1 p primary stream at station 1 
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Thrust Augmenting Ejector 



Fig. 1 E7/Ejector configured aircraft. 



Fig. 2 Constant area mixing region. 
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Element k 
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TOTAL EJECTOR THRUST (Newtons) 



Inner jet boundary 

Fig. 5 Jet boundary streamline. 



TIME (milliseconds) 


Fig. 6 Transient flow thrust prediction for various C i values. 
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